Quark Coalescence for Charmed Mesons in Ultrarelativistic Heavy-Ion Collisions 
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We investigate effects of charm-quark interactions in a Quark-Gluon Plasma on tfie production of 
D and J/tp mesons in high-energy heavy-ion collisions. Employing a previously constructed coales- 
cence model that successfully reproduces the transverse momentum (pt) spectra and elliptic flow 
{v2{pt)) of light hadrons at RHIC from underlying light-quark distributions at the phase transition 
temperature Tc, D-meson and J/ip pT-spectra are evaluated. For the charm-quark distributions, 
we consider two limiting scenarios: (i) no rescattering, corresponding to perturbative QCD (pQCD) 
spectra and (ii) complete thermalization including transverse expansion. We find that D-mesons 
acquire a minimal 112 inherited from their light-quark content and corresponding semileptonic de- 
cay spectra of single electrons practically preserve the «2 of the parent particles, exhibiting marked 
differences between the pQCD and thermal scenarios for pr > 1 GeV. Likewise, the pT-spectra and 
yields of J/ip's differ appreciably in the two scenarios. 

PACS numbers: 25.75.-q,25.75.Dw,25.75.Ld 



I. INTRODUCTION 

Collisions of heavy nuclei at the Relativistic Heavy-Ion 
Collider (RHIC) are expected to provide conditions favor- 
able for the creation of a deconfined and chirally symmet- 
ric Quark-Gluon Plasma (QGP). Among the promising 
probes of this phase are hadrons containing charm (c) 
quarks. On the one hand, the abundance of observed 
J/i/j mesons was predicted to be sensitive to QGP for- 
mation P] due to Debye screening of the heavy-quark 
potential, leading to suppression of its production. On 
the other hand, spectra of open-charm states (mostly D- 
mesons) are believed to encode valuable information on 
charm-quark reinteractions in the hot and dense medium. 
Due to the relatively large charm-quark mass, rUc ~ 1.5 
GeV, cc production ispresumably dominated by primor- 
dial N-N collisions |2| , whereas thermalization of their 
momentum distributions is still an open issue, with im- 
portant ramifications for charmonium production. 

At RHIC energies, with an expected 10-20 cc pairs per 
central A-A collision, statistical models predict a sub- 
stantial regeneration of J/^p^s by recombination of c and 
c quarks close to the phase transition J., J,, ^ S 13 ■ 
These estimates reside on the assumption that charm 
quarks are in thermal equilibrium with the surrounding 
medium. Recombination approaches have also been ap- 
plied before with fair phenomenological success in study- 
ing flavor dependencies of open-charm production in el- 
ementary p-N and tt-A^ collisions 0, 0, 0|. Again, 
the heavy quark serves as a probe for the (nonther- 
mal) hadronization environment via coalescence-type 
processes with surrounding valence or sea quarks. Fi- 
nally, parton coalescence within a hadronizing QGP has 
recently been put forward as a mechanism for light 
hadron production at intermediate pt — 2-6 GeV in 
heavy- ion collisions at RHIC. The observed "anomalous" 
p/tt ratio of ~ 1, as well as the apparent "constituent- 
quark" scaling of the hadron V2, are naturally accounted 
for within this framework [H [ll El, El El El . 



In this paper, we employ quark coalescence to evaluate 
spectra of open- and hidden-charm mesons at RHIC, with 
the objective to address the following issues: (a) What is 
the sensitivity of the J/^ abundance and px spectrum to 
the underlying momentum distribution of charm quarks? 
(b) How does the interplay of the charm- and light-quark 
distributions translate into the pT-spectrum and ellip- 
tic flow of Z3-mesons? Point (a) lifts the assumption of 
complete c-quark thermalization common to statistical 
models, which is important to discriminate regeneration 
from suppression mechanisms. A flrst study of this kind 
has been performed in Ref. where c-quark distribu- 
tions from pQCD have been implemented into a kinetic 
rate equation solved in the background of an evolving 
QGP. For c-quark rapidity densities of dNcc/dy 2.5 in 
central Au+Au collisions at RHIC, the final J/ip number 
has been found to deviate from scenarios with thermal- 
ized c-quarks by ^ 50%. Our results to be discussed 
below differ from these estimates. Concerning point 
(b), the recent study of Ref. IS] has shown that single- 
electron Pt spectra from decays of PYTHIA-generated 
and hydrodynamic D-meson distributions cannot by dis- 
criminated by current PHENIX data 19], extending to 
~ 3 GeV. In the present work, Z?-mesons are formed 
from c-quark coalescence with thermal light quarks at 
hadronization, where the latter also carry the collective 
expansion characteristics (radial and elliptic flow) that 
underlies a satisfactory description of pion and baryon 
spectra in Refs. [13 . Iisl ]. An important point in the 
present study will be the evaluation of the elliptic flow of 
both ZJ-mesons and their electronic decay spectra. For 
both (a) and (b) we will consider two limiting cases for 
the c-quark momentum spectra, i.e., (i) pQCD distribu- 
tions from PYTHIA, representing no rescattering, and, 
(ii) complete thermalization including collective expan- 
sion. 
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II. COALESCENCE INTO CHARMONIUM 

Let us start by recalling the basic elements of the coa- 
lescence model. In the Wigner formalism, the pertinent 
PT'-spectrum of a meson M takes the form 
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with d<7i denoting space-like hypersurface elements. 
fq,q{xi,pi) are invariant distribution functions of anti- 
/quarks depending on their space-time position and 4- 
momentum, and including spin and color degeneracy. 
The statistical factor accounts for the probability of 
forming a colorless meson of given spin from the under- 
lying quark color and spin (e.g., 3^=1/12), whereas 
encodes the dynamical part of the process. In the follow- 
ing, we will adopt plane-wave single-quark distribution 
functions, and neglect binding-energy corrections to /a/, 
in which case it becomes a Wigner distribution function. 
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FIG. 1: Transverse-momentum spectra of charm quarks in 
central 200 AGeV Au-|-Au collisions at midrapidity for dif- 
ferent scenarios: static thermal (dash-dotted line), thermal 
plus transverse flow (solid line) and primordial pQCD (dashed 
line). 

Since the measured hadron distributions in heavy-ion 
collisions at RHIC energies are approximately uniform 
around midrapidity, we will assume this also for the anti- 
/quark distribution functions, fq,q, with an additional 
Bjorken-type correlation between longitudinal coordinate 
and momentum ^0] . In the transverse direction, we fol- 
low Ref. [l^ by employing light-quark momentum distri- 
butions that are a combination of thermal spectra with 
radial flow and quenched minijets according to the GLV 
approach 21] with an opacity parameter L/X = 3.5. For 
the transverse momentum distribution of charm quarks. 



two possibilities will be investigated: (i) pQCD spectra 
as obtained from the PYTHIA event generator |23| and, 
(ii) thermal spectra at a temperature of T=170 MeV with 
radial flow. For the latter case, we take as an upper limit 
the same flow profile as determined in Ref. 12] for light 
quarks, i.e., a linear increase with the position in the 
transverse plane. 



/3(r) = /3n.ax;^ 



(2) 



with i?— 8.1 fm, the radius of the firecylinder at 
hadronization, and /3max=0.5 or 0.65, representing in- 
eherent uncertainties within our approach. Since heavy 
quarks are expected to suffer less energy loss in the QGP 
than light quarks jU0], we neglect charm-quark energy 
loss in the PYTHIA spectra in present study. In Fig. ^ 
we show by the dashed line the charm-quark p^-spectrum 
from the PYTHIA event generator and by the solid line 
that from the thermal scenario with collective flow. For 
comparison. Fig. also contains a thermal c-quark spec- 
trum without collective flow (dash-dotted line). 

The wave functions of the charmonium states are taken 
to be of Gaussian form in coordinate space, implying that 
the corresponding Wigner function is a Gaussian as well, 
i.e., 



/m = 8 exp{-x^/a^) exp(-gV2) 



(3) 



with X = X1—X2, and = (mi —7712)^ — {pi —P2Y, where 
X and p are 4-coordinate and -momentum (see Ref. [T3 ] 
for details). The width of the Gaussians is related to 
the mean-square-radius r^^^=< > via r^j^g = 3(jV8. 
Typical values from Cornell- type potential models j25| 
amount to ri-ms=0.47 fm for J/ijj, 0.74 fm for Xc and 
0.96 fm for tp', which will constitute our baseline values. 

To illustrate basic features of the coalescence approach 
in a transparent way we first employ thermal c-quark 
spectra without radial flow. Under these conditions, the 
momentum spectrum, Eq. can be evaluated analyt- 
ically upon invoking the non-relativistic Boltzmann ap- 
proximation for the c-quark distribution functions, suit- 
able for charmonia at low and moderate px- One finds 
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(4) 

with Ncc the total number of cc pairs within the rapid- 
ity range \y\ < 0.5. The bulk system is characterized 
by the proper time r, volume Vc and temperature Tc of 
the fireball at hadronization. Typical values for these 
quantities are r=4.3 fin^c=170 MeV, and 14=900 fm^, 
as determined in Ref. [l2l| from numerical evaluations in 
the light-quark sector, which, after inclusion of trans- 
verse flow, reproduce well pT-spectra of pions, protons, 
and kaons in central Au-|-Au collisions at 200 AGeV. The 
main features of Eq. Q are: (a) a weak dependence 
on the hadronization time as long as r > 2-3 fm/c and 
rims < 0.5 fm; (b) a ratio o-^/K ~ ^hadron/Vc of hadron 
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eigenvolume over fireball volume, characteristic for co- 
alescence (note, however, that the entailing increase of 
J 1^ production with a is an artifact related to the ne- 
glection of binding energy effects; we will return to this 
issue below); (c) a dependence on the number of par- 
ticipants, A, as N^/Ncc oc Ncs/V^ oc A'^^^/A oc A^/^, 
consistent with statistical models 0, 0, IE ■ 



TABLE I: Number of charmonia produced from charm-quark 
coalescence at mid-rapidity in central collisions of Au+Au at 
^ = 200 GeV. 









Thermal 2.6x10"^ 


0.2x10"^ 


0.7x10-* 



PYTHIA 0.8x10^^ 0.5x10^* 0.2x10^'' 



III. CHARMONIUM SPECTRA 



The numerical evaluation of Eq. |^ is performed along 
the lines of R.ef . [T^ . For each scenario, the total number 
of cc pairs is fixed to Ncs = 2.5 over one unit of rapidity. 
This number follows from an extrapolation of charm pro- 
duction at fixed target energies to central Au-|-Au colli- 
sions at 200 AGeV, corresponding to leading-order pQCD 
calculations upscaled by an empirical X-factor of ^--^5 26j, 
as well as the number of N-N collisions. With a charm- 
quark mass of TTLc = 1.4 GeV, the resulting open-charm 
cross sections are consistent with first indirect measure- 
ments via single electrons in 130 and 200 AGeV Au-|-Au 
by PHENIX jl^ at all centralities. 




FIG. 2: J/tp transverse-mass spectra for central Au-|-Au col- 
lisions at ^/s — 200 AGeV from charm-quark coalescence: 
circles correspond to PYTHIA c-quark distributions, while 
filled and open triangles follow from thermal spectra with ra- 
dial flow velocities of /3max=0.5 and 0.65, respectively. Lines 
represent exponential fits where appropriate. 

Fig. 121 displays our results for J/^ mT-distributions 
(without feeddown corrections); since, due to its large 
mass, the J/tp is more sensitive to variations in the col- 
lective expansion then light hadrons, we indicated in- 
herent uncertainties by using two values for the radial 
fireball-surface velocity at Tc, /3max=0.5 and 0.65. We 
find that both pQCD-based (circles) and thermal c-quark 
spectra with /3inax=0.5 (filled triangles) exhibit approx- 
imately uniform slopes, which are quantified by the fit- 



ted fines to be Tsiope=720 MeV and Tsiopc=350 MeV, 
respectively (for /3jiiax=0.65, Tsiopc=550 MeV in the re- 
gion rriT — Tajj^ = — 1.5 GeV). The average < pr > 
decreases from 1.95 GeV for the pQCD-based case to a 
range of 1.29-1.60 GeV for the thermal one. Further- 
more, the total Nj/^ obtained by employing c and c 
spectra from pQCD is smaller than that from the ther- 
mal model by about a factor of 3, as the flatter pQCD 
spectrum entails a substantial depiction of the probabil- 
ity in momentum space to form charmonia at low and 
moderate px where their yield is concentrated. Contrary 
to the single-electron spectra from heavy-meson decays, 
which appear to be relatively insensitive to the difference 
between pQCD and thermal-f flow open-charm distribu- 
tions ^3 (at least up to pT — 3 GeV, and including "con- 
taminations" from _B-meson decays), J /ip spectra exhibit 
a more pronounced sensitivity, which could enable it a 
rather direct window for studying the charm-quark be- 
havior in the QGP. 

Let us briefly address the question of feeddown cor- 
rections from Xc and ip' resonances. A rigorous treat- 
ment of their formation requires the inclusion of in- 
medium effects on binding energies, charm-quark masses 
and radii in the coalescence probability, which is a much 
more involved problem. First studies in this direc- 
tion, including formation and dissociation processes, have 
been undertaken in Refs. (t], H, ll^l- For simplicity, we 
take as a guideline thermal mass weights according to 
{niR/mj/^)^/^ exp{—{mR — mj/^,)/T) for their produc- 
tion relative to that for J ftp. The J/ip numbers corre- 
sponding to the spectra in Fig.|21are collected in Tabled 
together with estimated feeddown corrections. Despite 
the significantly larger radii of Xc and -0' states, their 
combined contribution to J/ip is still only ~10%. A more 
precise assessment is beyond the scope of this work. How- 
ever, we expect the relative suppression of the J/tp yield 
and its harder slope, as obtained with the pQCD charm- 
quark spectrum compared to the thcrmal+flow case, to 
be a rather robust result. It could prove valuable in, e.g., 
disentangling a transition from coalescence of soft c-c to 
hard ones (or even unsuppressed primordial charmonia) 
in pt, as well as in probing the degree of thermalization 
of charm quarks in the QGP. 
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IV. 73-MESON SPECTRA 

The standard coalescence approach becomes unrehable 
if the phase space density of either one (or both) of the 
quarks is no longer small, such as in bulk hadron produc- 
tion at low pt- For the case of interest here, i.e., D-meson 
spectra, we expect the following results to be trustwor- 
thy for Pt > 1 GeV. Again, our calculations proceed 
along the lines of Ref. |0| , with the same fireball param- 
eters as used above, i.e., a constituent light-quark mass of 
300 MeV, and the two scenarios for charm quarks (pQCD 
vs. thermal-|-flow) . The width of the Gaussian Wigner 
function for ZJ's is taken to render their radius to 0.6 fm. 
We have again weighted the production of D*'s with 
a thermal factor (tod* /mu)'^/^ exp(— (tod* — m£))/T) 
relative to that of D mesons. Semileptonic decays of 
ZJ-mesons have been computed assuming predominance 
of 3-body decays, Kve and K*ve^ with a phase-space 
weighting according to the weak matrix element. 

In Figl^l the upper curves display the coalescence re- 
sults for the pT-spectra of charmed mesons D, D*, Dg 
and D*. Also for Z3-mesons, the pQCD c-quark distribu- 
tions induce harder spectra (dashed line) than the ther- 
mal-|-flow ones, although somewhat less pronounced than 
for the charmonium case; the shaded bands in Fig|21cover 
the range of flow velocities from /3max=0.5 to 0.65). The 
lower curves show single- spectra from decays of all D- 
and i^s-mesons (including feeddown from D*'s). Similar 
to Ref. ^3 1 fi'^d that, on the basis of the current 
experimental accuracy, also within a coalescence frame- 
work, underlying pQCD c-quark spectra cannot really 
be distinguished from thermalized ones with flow. Our 
coalescence spectrum from thermal c-quarks with flow in- 
deed closel y re sembles the hydrodynamic Z?-meson spec- 
tra of Ref. fl3| (implying that, modulo charmed baryons, 
essentially all c-quarks recombine) . This is also true when 
comparing the pQCD-based c-quark coalesence with the 
fragmentation spectrum in Ref. above pr— 2 GeV. 
Below Pt — "2 GeV, however, there is a significant lack in 
the yield of the pQCD coalescence spectrum as compared 
to inclusive PYTHIA D-mesons, since no fragmentation 
contribution has been accounted for in the present study. 

Concerning the coalescence yields for Dg mesons (in- 
cluding feeddown from D*), we find for the two cases of 
thermal and pQCD c-quark spectra Dg/D ratios of 29% 
and 25%, respectively, with spectral shapes very similar 
to the inclusive ones of Fig. O These values are quite 
reminiscent to the 23% found within the coalescence ap- 
proach of Ref. 01 for p-p collisions at RHIC. 



V. CHARM AND SINGLE-ELECTRON 
ELLIPTIC FLOW 

Finally, we turn to the elliptic flow, t>2, of charmed 
particles, which is a well-established probe of the degree 
of thermalization in the (early) QGP as it vanishes in 
the absence of any rescattering. The coalescence model 




Pt (GeV) 

FIG. 3: Coalescence results for transverse-momentum spec- 
tra of all _D-mesons (upper curves) and pertinent electron de- 
cays (lower curves) in central 200 AGeV Au-|-Au collisions at 
midrapidity. Dashed lines are obtained from pQCD c-quark 
distributions, while shaded bands correspond to thermal spec- 
tra with flow velocities between 0.5 and 0.65. "Non-photonic" 
single-e* data (squares with both statistical and systematic 
errors) are from PHENIX O. 

has been shown to quantitatively describe [TsL ll^ |2^ 
the elliptic flow of light hadrons such as pions, protons, 
kaons and lambdas. Therefore, we believe that its perti- 
nent predictions especially for ZJ-mesons, and to a lesser 
extend for charmonia, are a rather reliable measure of 
charm-quark reinter actions. Even if c-quarks do not ex- 
perience any flnal-state interactions, a lower limit for the 
V2 of Z3-mesons from coalescence arises due to the V2 of 
their light-quark component. On the contrary, an up- 
per limit can be expected from the assumption that the 
charm quark has the same v^'^^ as the light quarks. 

Some specific features of charmed hadron V2(j>t) have 
already been pointed out within a schematic calculation 
in Ref. . In particular, it has been shown that even if c 
quarks acquire the same V2{pt) as light quarks, the large 
difference in the quark masses entails a smaller value of 
V2 for D-mesons than for J/ij] at same pt- This is so 
because the £)-meson momentum is largely determined 
by the c-quark, and an equal velocity of c and light quark, 
as required for coalescence, implies the light quark to be 
at low momentum where its V2 is relatively small. Here, 
we quantify this study by including more realistic quark 
distribution functions, especially through the inclusion of 
radial flow, which by itself affects light and heavy quarks 
in a different way, and extend it to the semileptonic decay 
electrons. 

Fig. 0] summarizes our elliptic flow results. For the up- 
per panel, we recall the underlying light-quark V2 (solid 
circles), which has been fit (solid line) to reproduce the 
experimental data for pions and charged hadrons (filled 
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FIG. 4: Elliptic flow in minimum bias Au+Au collisions at 
\/s= 200 AGeV. Upper panel: input parameterizations for 
light (solid circles) and charm quarks (open circles) together 
with the «2 for pions (solid line) from coalescence plus frag- 
mentation, as well as experimental data from PHENIX for 
(filled sauaresl|30l| and charged particles (open squares) [sH - 
Lower panel: D-mesons (solid lines, lower: using pQCD 
charm-quark distributions, upper: for thermal-l-flow c-quark 
spectra at T^) and respective decay electrons (diamonds), as 
well as J/ip mesons (dash-dotted line). 

and open squares, respectively) in Ref. together with 
a charm-quark V2 (open circles) saturating at the same 
maximal value (and with equal radial flow). The lower 
panel contains the ensuing results for D-mesons (solid 
lines) and their decay electrons (diamonds), as well as 
J/ip^s (dash-dotted line). For Z?- mesons, the elliptic flow 
from the two scenarios of complete thermalization and 
no reinteractions of the c-quarks deviates by about a 
factor of 2 for pt>1-5 GeV. The decay are found 
to essentially preserve the V2 of their parent particles, 
which makes it a highly interesting observable for up- 
coming measurements. We note that B-meson contribu- 
tions have not been included in present study. If these 
become signiflcant for momenta above 2 GeV, cau- 
tion needs to be exercised in drawing conclusions about 
c-quark elliptic flow. 

For the elliptic flow of J/V''s, which entirely stems from 
(thermalized) c-quarks, we observe a significant delay of 
its increase with pT as compared to that for Z?-mesons, 
caused by the radial flow. Indeed, if the latter is ab- 



sent, the V2 of J/ip's saturates faster than for D-mesons, 
see, e.g., Ref. [23 ■ We furthermore recall that, on the 
one hand, the V2 of J/V''s is zero for coalescence between 
pQCD charm quarks. On the other hand, we have not 
taken into account any primordially produced J/i/^'s, for 
which suppression mechanisms, if active early in the evo- 
lution, can induce a nonzero V2 113 ■ However, the mag- 
nitude of such effects (which can arise from both nuclear 
absorption and QGP dissociation) has been estimated to 
be rather small, up to ~2% (compared to 10% or more 
for coalescence between thermalized c-quarks). Never- 
theless, in a more complete description, one expects an 
increasing fraction of (unsuppressed) primordial charmo- 
nia towards higher px, thus reducing V2- 



VI. SUMMARY 

Within a coalescence approach as successfully ap- 
plied earlier in the light-quark sector, we have evaluated 
transverse-momentum dependencies of charmed hadrons 
in central heavy- ion reactions at RHIC. For the charm- 
quark distributions at hadronization we have considered 
two limiting scenarios, i.e., no reinteractions (using spec- 
tra from PYTHIA) and complete thermalization with 
transverse flow of the bulk matter. The resulting J/ip 
{rriT-) spectra differ in slope by up to a factor of 2 (harder 
for pQCD c-quarks), and the integrated yield is about a 
factor of 3 larger in the thermal case. For I?-mesons, 
we found that the difference in the slope parameters of 
the pT-spectra in the two scenarios is less pronounced, 
but their elliptic flow is about a factor of 2 larger for 
Pt>1-5 GeV in the thermalized case. The elliptic flow 
pattern of _D-mesons was found to be essentially pre- 
served in the single-electron decay spectra, rendering the 
latter a very promising observable to address the strength 
of charm reinteractions in the QGP. The present study 
can be straightforwardly generalized to charmed baryons 
(Ac), which may serve as a complimentary probe for 
charm-quark reinteractions in the QGP. 
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